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ABSTRACT 


This  paper  develops  a  non-linear  dynamic  capacitated  network  model  for 
planning  the  movements  of  pilgrims  in  the  Hajj,  one  of  the  world's  largest 
mass  movements,  according  to  religious  ritual,  which  would  assist  in  minimizing 
traffic  congestion  and  the  overcrowding  of  the  holy  sites.  A  new  non-linear 
representation  of  congestion  with  convenient  mathematical  properties  is  made. 
The  model  is  effective  in  producing  quantitative  and  qualitative  background  for 
general  policy  decisions  on  the  Hajj  transportation. 
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A  NONLINEAR  CONGESTION  NETWORK  MODEL  FOR  PLANNING 
INTERNAL  MOVEMENTS  IN  THE  HAJJ 


by 

A.  Charnes,  S.  Duffuaa,  and  A.  Yafi 


1 . 0  Introduction 

The  purpose  of  this  paper  is  to  identify  issues  in  the  Hajj  internal 
transportation  network  and  model  the  situation  in  order  to  draw  general  policy 
decisions  in  management  of  this  unique  situation.  The  most  important  constraints 
in  this  case  are  the  religious  ones,  since  they  define  the  transportation 
and  the  movement  of  the  people  from  one  place  to  another.  Any  model  must  reflect 
these  constraints  in  order  to  be  implementable.  These  religious  constraints  are 
represented  as  time  and  space  constraints.  A  temporal  capacitated  network  is 
developed  to  represent  the  situation,  with  the  objective  to  minimize  congestion 
of  traffic  on  roads  and  overcrowding  of  the  holy  places.  This  objective  is 
chosen  to  enable  the  visitors  to  perform  their  duties  as  easily  as  possible. 

In  section  2  we  describe  the  Hajj  situation  briefly,  and  in  section 
3  we  identify  issues  to  be  considered  in  the  modeling  process.  In  section  4 
we  develop  a  linear  network  model,  and  in  section  5  we  present  a  numerical 
example  for  the  Hajj  situation.  In  section  6  we  extend  the  model  to  a  non-linear 

model  to  properly  account  for  congestion  via  a  non-linear  function  and  we 
draw  conclusions. 

2.0  Description  of  the  Hajj 

The  Hajj  is  an  annual  meeting  for  about  two  million  Muslims  to  achieve 
one  of  the  five  pillars  of  Islam.  The  Hajj  involves  several  visitations  to 
several  holy  sites  and  it  occurs  every  354  days,  since  it  is  dated  by  the 
Islamic  Lunar  year.  The  Hajj  process  in  this  paper  will  be  divided  into  three 


phases.  Phase  one  mainly  consists  of  air  and  sea  transportation  of  foreign 
population  to  Jeddah,  and  also  inland  transportation  of  foreign  and  native 
population.  The  duration  of  this  phase  is  up  to  the  eighth  of  the  last  month  of 

the  Lunar  year,  by  that  time  all  the  pilgrims  have  gathered  in  one  of  the 
following  towns:  Jeddah,  Makkah,  Muna  or  Medina. 

Phase  two  starts  from  the  morning  of  the  eighth  of  the  last  month 
in  the  Lunar  year  up  to  the  end  of  the  Hajj  and  the  arrival  of  the  pilgrims 
at  the  ports  to  leave  for  home.  This  phase  includes  part  of  the  in-land 
transportation  from  Jeddah  and  Medina  to  Makkah  and  all  the  movement  between 
the  holy  sites  up  to  the  completion  of  the  Hajj.  Its  duration  is  roughly  up 
to  the  nineteenth  or  the  twentieth  of  the  month.  Phase  three  consists  partly 
of  the  transportation  between  Makkah  and  the  ports.  It  also  includes  the 
departure  of  all  prilgrims  for  home.  In  this  paper  we  deal  with  phase  two, 
the  "internal  transportation  network  of  the  intramovement  of  the  pilgrims", 
which  constitutes  the  logistics  of  the  religious  process. 

2.1  Phase  Two  of  the  Hajj 

In  general,  the  Hajj  season  begins  at  the  start  of  the  tenth  month 
of  the  Muslim  lunar  calendar  year  and  ends  about  the  twentieth  of  the  Dihu'L- 
Hijjah  (the  twelfth  and  final  month  of  the  lunar  year).  Most  of  the  pilgrims 
arrive  in  the  last  fifteen  days  of  the  season.  According  the  our  model, 
phase  two  starts  the  morning  of  the  eighth  day  of  the  last  month  of  the  lunar 
calendar  year.  By  this  time  all  pilgrims  have  already  arrived  in  Saudi  Arabia 
and  they  are  either  at  Makkah,  Muna,  Jeddah  or  Medina  (see  dagram  in  Figure  I). 

It  is  a  ritual  that  before  arriving  at  Makkah  every  pilgrim  has  to 
wear  the  garment  of  iharm  (restriction) ,  is  forbidden  to  hunt,  argue,  cut  his 
hair,  clip  his  nails  or  engage  in  any  sexual  activities.  Also  on  arrival  at 


Makkah  (45  miles  east  of  Jeddah)  each  pilgrim  has  to  make  the  greeting  tawaf, 
the  prescribed  seven  counterclockwise  circumambulations  of  the  Kaaba  (a  black 
room  in  the  center  of  the  holy  mosque  built  by  the  prophet  Ibraham).  With 
the  tawaf  they  perform  the  Sa'y--making  seven  trips  between  the  hills  of 
Safa  nad  Marwah.  The  path  of  Sa'y  is  enclosed  in  a  long  gallery  which  is 
part  of  the  holy  mosque. 

Also  on  the  eighth  day,  those  pilgrims  who  are  not  already  in  Muna, 
four  miles  east  of  Makkah,  move  through  Makkah  to  Muan  for  the  essential 
final  days  of  the  Hajj.  The  next  day  (on  the  ninth)  everybody  has  to  be  in 
Arafat,  eight  miles  east  of  Muna,  to  perform  the  Standing,  the  central  ritual 
of  the  Hajj.  The  duration  of  the  Standing  is  from  noon  of  the  ninth  until 
sunset. 

On  the  eve  of  the  tenth,  the  pilgrims  leave  Arafat  and  stop  at 
Muzdalifa,  a  place  between  Muna  and  Arafat.  At  Muzdalifa  they  collect 
pebbles  to  throw  at  the  three  "Satan's  stoning  points"  in  Muna  during  the 
following  days.  These  points  symbolize  the  force  of  evil.  Those  who  leave 
for  Makkah  start  doing  one  of  the  pillars  of  the  Hajj  Tawaf  El  ifada,  the 
post-Arafat  tawaf  done  in  the  same  manner  as  the  greeting  tawaf.  After 
finishing  this  tawaf  pilgrims  can  put  off  iharm  (restriction)  and  return  to 
Muna  to  finish  the  Stoning  of  the  devil.  After  finishing  the  Stoning,  the 
Hajj  is  complete. 

All  pilgrims  after  the  Hajj  perform  a  farewell  tawaf  and  some  of 
them  leave  for  Medina  to  visit  the  prophet's  grave.  Visiting  Medina  can  be 
done  before  the  Hajj  starts  and  those  who  did  it  before  the  Hajj  leave  via 
Jeddah  or  by  inland  routes  for  home.  Figure  1  shows  all  the  holy  sites  and 
the  routes  of  the  pilgrims'  movements. 


Figure  1 


A  Diagram  for  the  Islamic  Holy  Sites  and  the  Flow  of  Pilgrims 


From  Asia, 
Jordan,  Syria 


countries  From  Yemen 


4 

Arafat 


>  The  solid  lines  represent  the  internal  movement 

->  The  broken  lines  represent  the  external  movement 
for  coming  to  Hajj  and  leaving 


3.0  Issues  to  be  Considered  In  the  Model 

The  most  important  elements  here  are  the  timing  and  place  of  each 
of  the  holy  practices.  This  determines  the  flow  of  people  from  one  place  to 
another.  We  notice  from  section  two  that  there  is  some  flexibility  in  the 
timing  of  most  of  the  pillars  of  the  Hajj.  The  model  must  be  able  to  take 
full  advantage  of  that  in  order  to  make  the  Hajj  as  smooth  as  possible. 

When  the  pilgrims  travel  to  Makkah  they  are  in  large  numbers  and 
there  is  overcrowding  of  the  roads  which  lead  to  Jakkah.  On  arrival  at 
Makkah  they  perform  the  greeting  tawaf.  There  they  need  some  scheduling  to 
ease  overcrowding  at  the  Kaaba. 

The  central  ritual  of  the  Hajj  is  the  Standing  at  Arafat.  This 
is  a  large  place  and  can  take  all  the  pilgrims,  but  on  the  way  to  and  from 
Arafat  the  big  crunch  occurs.  Then  the  people  shuttle  between  Makkah  and 
Muna  until  they  finish  their  Hajj. 

In  scheduling  events  in  the  Hajj,  the  following  must  be  noted. 
First,  the  greeting  tawaf  can  be  done  any  time  before  the  morning  of  the 
ninth  of  the  Hajj  month.  The  post-Arafat  tawaf  can  be  performed  at  any  time 
from  the  eve  of  the  tenth  to  the  evening  of  the  twelfth  of  the  Hajj  month, 
but  most  people  like  to  do  it  earlier  to  finish  their  Hajj.  Secondly,  the 
Stoning  of  the  Satan  is  flexible  within  two  or  three  days,  i.e.,  from  the 
tenth  up  to  the  twelfth  of  the  Hajj  month.  Third,  the  Sacrifice  of  the 
sheep  can  be  done  any  time  from  the  morning  of  the  tenth  until  the  evening 
of  the  twelfth,  but  most  people  do  it  on  the  tenth.  Fourth,  most  people 
catch  the  'Id  prayer  at  the  Sacred  Mosque  on  the  morning  of  the  tenth. 

Fifth,  the  farewell  tawaf  can  be  done  any  time  before  leaving  Makkah.  Sixth, 
the  visitation  to  Medina  can  be  done  before  or  after  the  Hajj  and  there  are 
no  time  constraints  on  it. 


The  model  should  take  advantage  of  the  flexibility  in  the  timing 
of  the  Hajj  rituals,  but  also  should  allow  room  for  pilgrims'  preferences  in 
choosing  the  time  to  do  rituals  if  the  models'  specifications  are  to  have  any 
chance  of  being  implemented. 

4.0  Statement  of  the  Model 
4.1.  The  Constraints  Set 

The  following  constraints  set  of  the  Hajj  internal  movement 
is  shown  in  Figure  2.  Figure  2  represents  movements  of  pilgrims  by  half-days 
from  the  eighth  of  the  Hajj  month  till  the  nineteenth.  Each  half  a  day  is 
considered  as  a  period,  so  the  model  consists  of  twenty-four  periods.  We 
think  half  a  day  is  a  reasonable  period  for  general  policy  decisions,  but  for 
close  monitoring  of  the  pilgrims'  movements  the  model  can  be  represented  in 
terms  of  shorter  periods,  i.e.,  in  terms  of  one  hour  period.  This  will  expand 
the  model,  but  since  the  model  is  of  network  type,  it  can  be  solved  easily. 

Figure  2  gives  a  complete  representation  of  the  Hajj  internal 
movements.  There  are  two  sets  of  arcs  in  Figure  2,  one  set  represents  movement 
of  pilgrims  within  periods,  the  other  set  represents  pilgrims  who  are  staying 
in  the  same  town  for  the  next  period. 

The  upper  bounds  on  the  arcs  in  the  model  should  be  considered 
very  closely  and  the  following  considerations  would  be  very  useful: 

(1)  the  time  of  the  period  evening  or  morning, 

(2)  the  lodging  capacity  of  each  town, 

(3)  the  food  and  public  facilities  of  each  location, 

(4)  the  health  facilities, 

(5)  the  availability  of  buses,  cars,  between  any  two  towns  on 
each  route,  and 


(6)  the  timing  of  each  ritual. 


FIGURE  2 

NETWORK  REPRESENTATION  For  THE  HAJJ 
Internal  MOVEMENT 
Each  Half  A  day  is  A  period 
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s  represent  Jeddah,  2  h represent  Makkah.  3  =  represent  Muna,  4  =  represent  Arafat. 
5  s  represent  Medina.  3MZ  =  represent  Mozdaiifa 
0  -  Dav.  N  s  Night 


=  number  of  pilgrims  going  from  location  i  to  location  j  in 


period  t; 

y.jj(t,t+l)  =  number  of  pilgrims  going  from  location  i  to  location  j 
during  period  t; 

y^(t,t+l)  =  number  of  pilgrims  staying  at  location  i  during  period  t; 
[In  this  case  we  do  not  consider  in  detail  cross  travel  between  -iods 
from  one  site  to  another,  since  one-half  day  is  more  than  suffir  nt  travel 
time  between  any  two  sites;  therefore,  flexibility  in  the  timing  o'  ituals 
enables  the  model  to  accommodate  cross  travel.] 

K.  .  =  the  maximum  number  of  pilgrims  who  can  travel  easily  at 

*  J 

day  times  from  location  i  to  location  j,  i.e.,  this  is 
the  capacity  of  route  (i,j) 

L. -  =  the  maximum  number  of  pilgrims  who  can  travel  easily  at 

1  J 

night  time  on  route  (i,j) 

Ln-  =  lodging  capacity  of  town  i, 

K..  =  food  capacity  of  town  i, 

b..  =  the  minimum  of  the  public  facilities'  capacities  at  town  i. 

Let  1  be  Jeddah,  2  be  Makkah,  3  be  Muna,  3.5  be  Muzdalifa  ,  4  be  Arafat,  and 
5  be  Medina. 

P.j  =  number  of  pilgrims  at  town  i  at  the  beginning  of  the  first 

period. 

The  following  are  the  equations  defined  by  the  network  in  Figure  2.  The 
numbers  to  the  left  of  the  equations  represent  the  period  and  the  node,  i.e. 
1.5  represent  the  equation  for  period  one  at  node  5,  or  the  equation  for 


'w 


The  following  periods  19,  20  up  to  24  have  the  same  pattern  or  structure 


as  18  so  we  will  have  just  copies  of  the  equations  in  18.  The  following  are 
the  capacity  constraints  on  the  routes  and  the  holy  places. 

0  <  x.  ■  (t)  <  K. .  if  t  is  odd 

l  <J  •  J 

0  <  x..(t)  <  L..  if  t  is  even 

i  o  i  J 

0  <  Vi  j 

4.2.  Objective  Function 

One  of  the  objectives  to  be  considered  is  space  and  routes  utilization. 
The  objective  is  of  goal  type  to  ease  crowding  as  much  as  possible,  i.e. 

(4.2.1)  Max  |  min  (K^ ,L^ ,b - )  -  yi i (t,t+l)|  +  EE 1  Ku  ■  xij(t)l 


t  iJ 
odd 


+  EE 1  Li j  -  *,■  j 

t  ij 
even 


i(t)l  *EE|LiJ  - 

t  ij 


Maximizing  this  function  alone  will  allow  for  overcrowding  of  holy  sites  and 
we  will  be  maximizing  a  convex  function,  but  with  the  capacity  constraints 


xij(t>  *  Lij  •  xij(t)  4  Kij  and 
y^t.t+l)  <min(K.,  L . ,  b.)  ,  y...(t,t+l)  < 

that  will  not  happen  and  the  objective  function  will  be  transformed  into 
the  following  equivalent  linear  function. 

(4.2.2)  Mfn^E  yii(tlt+1)  +EE  xij(t)  +  EEE 

t  i  t  i  j  t  i  j  t 

Hence  the  model  for  the  Hajj  internal  movement  will  minimize  the  linear  function 

in  (4.2.2)  subject  to  all  the  constraints  in  section  4.1. 

The  above  model  could  be  updated  if  conditions  on  the  routes  and 
the  holy  sites  are  changed.  It  could  be  formulated  in  terms  of  shorter 

periods,  for  example  on  an  hourly  basis,  if  detailed  monitoring  of  the  Hajj 
internal  movement  is  needed. 


Another  objective  function  to  be  considered  is 
n 

Hinyyj  min(Ki,L.tb.)  -  y.(t,t+l)|  +EE|Kij'xu(t)l 
t  i=l  t  ij 

odd 

+EE  i  Lij  -  +  EE1  Lij  •  yu(t-ul)i 

t  ij  t  ij 

even 

This  is  another  type  of  a  goal  type  objective  and  this  with  the 
constraints  in  section  4.1  can  be  transformed  into  a  larger  pure  network. 

Other  objectives  and  constraints  may  also  need  consideration,  e.g., 
the  role  of  individual  contractors  ( "Mutwafeen")  may  need  more  explicit 
consideration. 


5.0  Examples  and  Extensions 

The  following  example  is  a  reasonable  representation  for  the  Hajj 
situation  in  the  year  1980.  The  example  was  constructed  using  [4],  a  study  done 
by  the  Hajj  Research  Center  at  Jeddah.  Using  this  study  and  some  personal 
judgement,  the  following  data  on  routes  and  town  capacities  were  specified. 


Routes 

Jeddah  to  Makkah  0 
Makkah  to  Muna  0 
Makkah  to  Medina  0 
Makkah  to  Arafat  0 
Muna  to  Arafat  0 
Muna  to  Makkah  0 
Medina  to  Makkah  0 


Day 

Night 

to 

200, 

,000 

0 

to 

150,000 

to 

500, 

,000 

0 

to 

400,000 

to 

200, 

,000 

0 

to 

150,000 

to 

900, 

,000 

0 

to 

700,000 

to 

800, 

,000 

0 

to 

700,000 

to 

250, 

,000 

0 

to 

200,000 

to 

100, 

,000* 

0 

to 

150,000 

*when  2-way  traffic 


2 


Capacity  Of  No.  of  Pilgrims 


Jeddah 

300,000 

Makkah 

1,000,000 

Muna 

600,000 

Arafat 

2,000,000 

Medina 

1,000,000 

Px  =  200,000  ;  P2  =  800,000  ;  P3  =  400,000  ;  P&  =  200,000  ;  P4  =  0 

The  above  example,  with  some  variations  to  allow  for  pilgrim  increase 
was  solved  in  order  to  evaluate  the  model  and  its  representation  of  the  real 
situation.  To  achieve  this  evaluation,  the  dual  variables  of  the  model  were 
utilized.  The  dual  variables  give  an  evaluation  of  the  model  response  to 
traffic  volumes  and  routes  capacities. 

In  all  three  examples  the  dual  evaluations  point  to  the  route  Mina 
to  Mecca  at  period  four  as  pre-eminent  in  relief  of  congestion  since  it  was 
the  most  negative  dual  variable.  An  increase  in  the  capacity  of  this  route 
by  one  unit  will  result  in  127  units  reduction  in  the  objective  function 
(which  increases  congestion)  since  the  value  of  its  dual  variable  is  -127. 

The  examples  further  show  that  major  concern  should  be  given  to  the  Hajj  in¬ 
ternal  movement  in  period  four  and  five.  The  complete  solution  of  this 
example  is  in  Table  1  in  the  Appendix. 

The  other  negative  dual  variables  occur  for  the  routes  Jeddah  to 
the  Sink  and  Medina  to  the  Sink.  The  most  negative  dual  variable  is  associated 
with  the  capacity  of  the  route  Jeddah  to  the~Sink  of  value  -14.  That  means 
an  improvement  in  Jeddah  port  capacity  of  one  unit  will  result  in  a  reduction 
of  14  units  in  the  objective  function.  Using  the  magnitude  of  the  dual 


variables,  the  ports  can  be  ranked  in  terms  of  their  importance  in  reducing 
traffic  congestion.  The  Jeddah  port  is  the  most  important  to  the  Hajj  in¬ 
ternal  movement  and  Medina  is  next. 


The  second  and  the  third  examples  give  directions  where  improvement 
should  be  made  in  case  the  number  of  pilgrims  increase.  Both  examples  show 
that  the  Jeddah-Mecca  route  must  be  improved  to  increase  its  capacity  to 
handle  the  anticipated  increase  in  pilgrims  in  the  future. 

These  examples  give  some  insight  into  decision  making  for  Hajj 
scheduling  but  the  simple  objective  function  does  not  reflect  very  well  all 

the  possible  effects  of  levels  of  congestion  on  different  routes.  Hence,  they 
give  scheduling  plans  which  are  insensitive  to  routes  and  to  changes  in 

capacity.  A  more  appropriate  measure  of  congestion,  which  will  result  in 
better  distributions  of  traffic  between  the  routes,  is  needed. 

Such  a  measure  of  congestion  will  be  presented  in  the  next  section. 
It  results  in  a  non-linear  objective  function  network.  We  solve  it  by  piece- 
wise  linearization  and  shall  compare  its  results  with  those  of  this  model. 


6.0  Non-Linear  Network  for  the  Hajj  Internal  Movement 

In  the  previous  section  we  examined  a  linear  model  for  the  Hajj  Inter¬ 
nal  Movement.  The  linear  objective  is  not  an  appropriate  measure  of  congestion 
since  it  does  not  correctly  reflect  the  fact  that  rate  of  congestion  increases 
with  increases  in  closeness  to  capacity  on  a  route  segment.  In  this  section 
we  introduce  a  non-linear  measure  of  congestion  for  each  route  which  depends  on 
the  capacity  of  the  route  segment  and  the  density  of  traffic  on  it.  In  the 
next  section  we  will  elaborate  on  the  properties  of  this  measure. 
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Let  R.j j ( x j{t) )  and  R  — (y^ j(t,t+l))  be  convex  non-linear  functions 
which  measure  the  congestion  on  arcs  x^(t)  and  y^ft.t+l).  The  model  in 
section  4.2  is  now  altered  (only  in  the  objective  function)  to  be 

(6-°-u  -  z  z  «„€*„««>  *  z  z 

t  ij  t  ij 

5 

subject  to  ^  e.jXijO:)  +  ^  ei J.(t)y...(t,t+1)  =  0 
J=1  j 

i  =  1,2,. ..,5  ;  t  =  1,2,. ..,24  ;  x.j(t),  y.(t,t+l)  >  0 

where  the  incidence  numbers  c..  are  1,  -1,  or  zero.  The  model  in  (6.0.1)  can 

1  J 

be  replaced  by  the  following  capacitated  linear  model  wherein  we  have  approx¬ 
imated  the  R..(x..(t))  and  R..(y.  .(t,t+l))  functions  by  two  linear  segments. 

•  J  '  J  I J  I  J 

In  the  next  section  we  shall  develop  a  method  for  choosing  the  two  line  seg¬ 
ments.  Employing  such  an  approximation  (6.0.1)  becomes 


(6.0.2) 


ZEZ  Z *EZZ  H 

t  i  j  k=l  t  i  j  k=l 


24  2  5  2 

subject  to  EZ  eijxijk(t)  +  zz,  j(t)yijk(t,t+l)  = 
j  k  j  k 


i  1 *2 , . . . ,5  ;  0  ^  xijk(t)  ^  Aijk  ;  0  ^  yijk(t)  ^  Ajjk 
The  dual  of  (6.0.2)  is 


(6.0.3)  24  5  24  5  5  2  24  5  5  2 

Max  ZZv11-0  -ZZZZ  'i'fjk(t>4ijk  -ZZZZ"ijk(t)Ajk 

t=l  i  t=l  i  j  k  t=l  i  j  k 


subject  to 


£e1J*1(t)  -  "ijk  4  Rijk  V  J-  l-  “here  “"restricted.  »1Jk>  »iJ|t  > 

In  an  optimal  solution  to  (6.0.2)  the  dual  variables  4>*.,  w*...  and 

1  1 J  K  1 J  K 

provide  evaluators  for  the  initial  points  at  which  changes  in  "resistance"  to 
traffic  volumes  occurs.  Hence  a  means  is  provided  for  evaluating  possible 
changes  associated  with  those  features  of  the  network  in  terms  of  their  effects 
on  the  resulting  traffic  volumes  at  the  times  associated  with  the  movements. 

To  develop  (6.0.2)  we  nedd  the  R.^,  A..  ^  for  the  model.  The  next 
section  will  elaborate  on  the  congestion  measure  we  have  chosen  and  on  how 
the  and  A..^  are  determined. 

6.1  A  Convex  Non-linear  Measure  of  Traffic  Congestion 
The  measure  of  Congestion  of  traffic  should: 

1)  increase  with  increase  in  traffic  flow, 

2)  have  an  increasing  rate  of  increase  with  traffic  flow  increase  (i.e.  be 
a  convex  function), 

3)  be  zero  or  very  close  to  zero  at  zero  flow, 

4)  approach  infinity  when  flow  approaches  the  capacity  of  the  route. 

One  such  measure  of  traffic  congestion  is 


(6.1.1)  f(x,.)  = 


Ui j  ’  X  ij 


where  is  the  capacity  of  the  route  movement  (i,j,).  To  apply  (6.0.2) 
we  need  to  find  and  A^  and  what  the  relationships  are  between  these  A^ 
and  R-jjk*  Thus  we  break  the  non-linear  function  into  two  line  segments. 

In  order  to  exhibit  further  relationships,  in  figure  6-1  we  exhibit  simul¬ 
taneously  two  examples. 


For  any  given  two  routes  (i,j)  and  (r,s),  they  have  the  measures  of 


congestion  Rfj(x^j)  and  Rrs(*rs)-  We  break  the  function  R^x^)  into  two 


line  segments.  The  first  line  segment  is  the  line  that  joins  the  two  points 


C 


0, _ 1  \  and 


1 


We  now  show  how  from  determination  of  the 


uij 


,u..,  U..-U.., 

^  Ijl  1J  Ijly 


breakpoint  for  one  route  (i,j),  we  can  simultaneously  determine  the  break¬ 
points  for  all  other  routes  (r,s).  Using  the  concept  that  congestion  will 
be  the  same  if  traffic  density  relative  to  route  capacity  is  the  same,  we 
consider 


(6.1.2) 


where  x  is  the  flow  in  route  (i,j)  and  y  is  the  traffic  flow  in  route  (r,s). 
Using  (6.1.2)  we  have 


Thus  the  breakpoints  are  determined  immediately  by  the  route  capacity  and 
so  are  the  slopes  for  all  the  routes. 


(6.1.5)  The  Rrs,  =  /  1_-  1  \  *  U 


u  -u  ,  u 

rs  rsl  rs 


rsl 


6.1.1  THEOREM 


If  R..  ^  is  the  slope  of  the  first  line  segment  for 


Rrsl=U*ij  '  R1jl 


for  all  (r,s). 


Proof: 


R~*(xrs)  = _ L 


rs 


Urs-Xrs 


Rrsl  iS  the  slope  0f  the  1ine  J°inin9  and /UrO  '-rr _ 1 


'rs  J 


\ 


rsl1 


U  -U  , 
rs  rsl 


'rsl 


Urs-Ursl  Urs 


T  U 


rsl 


U  -  U  +  U  , 
rs  rs  rsl 


1 


U  ,/u  -U  Au  U  /u  -  U  A 
rsl (  rs  rsl]  rs  rsl  rs  rsl] 


,  then 


To  determine  the  second  line  segment  slope,  let 


M  = 


where 


U .  •  -  x* . . 
U  iJ 


U..  '  x*ii  =  n  for  all  (i,j).  Then 

I  J  •  J 


1 


Rij2  Uij  "  Uijl 


X*ij  "  Uijl 


Computing  R..^,  for  all  functions  on  the  routes,  putting  them  and  their 
associated  variables  into  (6.0.2),  we  abtain  a  linear  network.  This  model 
has  been  applied  to  the  Hajj  Internal  Movement  using  the  same  data  as  in  the 
linear  case,  and  solved. 


6.2  Results  and  Conclusions 

The  non-linear  measure  of  traffic  congestion  is  applied  to  the 
same  example  previously  solved  for  the  linear  case.  All  the  necessary 
piece-wise  linearization  is  performed  as  described  in  section  6.1.  The 
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optimal  solution  for  the  non-linear  case  depicts  the  Hajj  situation  in 
more  realistic  terms  by  picking  up  the  required  movements  between  Muna  and 
Makkah  in  periods  5,  6,  and  7  without  imposing  lower  bounds  on  traffic 
movements  on  routes. 

Moreover  we  see  that  the  evaluation  given  by  the  dual  variables 
for  the  non-linear  case  give  a  clear  indication  of  routes  which  might 
ease  traffic  congestion. 

The  model  shows  that  the  time  interval  of  concern  for  the  Hajj 
Internal  Movement  is  from  period  four  to  period  seven.  To  comprehend  the 
ways  of  easing  the  traffic  congestion  there  we  need  to  concentrate  on  the 
Hajj  Internal  Movement  from  period  three  to  period  eight  and  use  the 
results  of  this  model  as  a  guide  for  traffic  improvements. 

The  model  evaluators  clearly  signal  the  following  routes  as 
the  key  routes  in  easing  the  traffic  congestion.  The  route  with  the  most 
negative  dual  variable  is  Muna  to  Makkah  in  period  four.  The  second  route 
is  Muna  to  Makkah  at  period  5,  followed  by  Muna  to  Muna  between  periods 
five  and  six. 

The  model  shows  that  the  pilgrims'  time  of  departure  depends 
strictly  on  the  capacity  of  the  ports.  Upgrading  the  capacities  of  the 
ports  will  help  in  smoothing  the  Hajj  process.  Another  important  way  to 
ease  the  traffic  congestion  is  to  divide  the  pilgrims  into  two  groups: 
one  to  visit  Medina  prior  to  the  Hajj,  the  other  after  the  Hajj.  Then 
the  group  which  visits  Medina  after  the  Hajj  departs  from  Medina.  This 
will  ease  the  traffic  after  the  Hajj. 

Complete  optimal  solutions  for  the  linear  and  the  non-linear 
cases  are  shown  in  Tables  1  and  2  in  the  appendix.  The  optimal  solution 


TABLE  1 
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Optimal  Solution  for  the  Linear  Example 
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